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1 Introduction 



The DIRAC collaboration ^ at CERN has measured the lifetime of pionium 
in its ground-state, and the preliminary result yields t„^io = [3.1^Q y(stat) ± 
l(syst)] • 10^ ^^s A lifetime measurement of pionium at the 10% level allows 
one to determine the S-wave tttt scattering length difference jap — a§| at 5% 
accuracy. The measurement can then be compared with theoretical predictions 
for the S-wave scattering lengths ||] and with the results coming from 

scattering experiments [|| . Particularly exciting is the fact that this enterprise 
subjects chiral perturbation theory to a very sensitive test Q. New measure- 
ments are proposed for CERN PS, J-PARC and GSI These experiments 
aim to measure the lifetime of tt+tt^ and tt^K^ atoms simultaneously. 
In order to extract the scattering lengths from such future precision measure- 
ments, the theoretical expressions for the energy shifts and decay widths of the 
TT+vr^ and n^K^ atoms must be known to a precision that matches the exper- 
imental accuracy. Nearly fifty years ago, Deser et al. ^ derived the leading 
order formulae for the decay width and the energy shift in pionic hydrogen. 
Similar relations exist for tt+tt^ and n^K^ atoms which decay due to 

the strong interactions into 27r° and tt^K^ , respectively. Theoretical investiga- 
tions on the spectrum and the decay of pionium have been performed beyond 
leading order in several settings. Potential scattering has been used Q 
as well as ficld-thcoretical methods |l^, |l^, ^ In particular, the 

lifetime of pionium was studied by the use of the Bethe-Salpeter equation |]l9j 
and in the framework of the quasipotential-constraint theory approach . The 
width of the tt+tt" atom has also been analyzed within a non-relativistic effec- 
tive field theory ||2l|, ^ , which was originally developed for bound states in 
QED by Caswell and Lepage |24|. The non-relativistic framework has proven 
to be a very efficient method to evaluate bound state characteristics. It was 
further applied to the ground-state of pionic hydrogen |2^, ^ ^ and very 
recently to the energy-levels and decay widths of kaonic hydrogen [^8| . Within 
the non-relativistic effective field theory the isospin symmetry breaking cor- 
rections to the Deser-type formulae can be evaluated systematically. In Refs. 

|2^ , ^ ^dj the lifetime of pionium was evaluated at next-to- leading order 
in the isospin breaking parameters a ~ 1/137 and (m„ — mdf' ■ 
We presented in Ref. the results for the S-wave decay widths and strong 

energy shifts of tt+tt^ and tt^K^ atoms at next-to-leading order in isospin 
symmetry breaking. Further, for the lifetime as well as for the first two energy- 
level shifts, a numerical analysis was carried out. The aim of this article is to 
provide the details that have been omitted in Ref. Chiral perturbation 

theory (ChPT) allows one to relate the result for the width of the ir^K^ atom 
to the isospin odd nK scattering lengths ^ while the strong energy shift is 
proportional to the sum of isospin even and odd scattering lengths a J + 0-0- 
The values for and Cq , used in the numerical evaluation of the widths and 
strong energy shifts, stem from the recent analysis of ■nK scattering from Roy 
and Steiner type equations |Q. Within ChPT, the ttK scattering lengths have 
been worked out at one-loop accuracy p^ , |3^ , p6[ , and very recently even the 
chiral expansion of the 'kK scattering amplitude at next-to- next-to-leading order 
became available . Particularly interesting is that the isospin even scattering 
lengths a J depends on the low-energy constant Lg |3^, and this coupling is 
related to the flavour dependence of the quark condensate . 
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The paper is organized as follows: The general features of tt+tt and ir^K^ 



atoms are described in Section 2 
proach is illustrated in Section 3 



The non-relativistic effective field-theory ap- 
by means of the 'k~ K'^ atom. The discussion 



includes the Hamiltonian, the master equation, and the matching to the rela- 
tivistic ttK amplitudes. In Sec. 0, we present the results for the decay widths 
and strong energy shifts of the n and 11^71^ atoms at next-to-leading order 
in isospin symmetry breaking. The role of transverse photons is discussed in 
Section |[ Transverse photons do contribute to the electromagnetic part of the 
energy shift. The pure QED contributions have been worked out a long time 
ago , based on the Bethe-Salpeter equation |}40|, the quasipotential approach 
[^H ^ and an improved Coulomb potential ||43| . We reproduced this result 
within the non-relativistic framework. We further estimate the contributions 
from transverse photons to the decay width of the atom and show that 

they vanish at next-to-leading order in isospin symmetry breaking. The con- 
tributions generated by the vacuum polarization of the electron |2^, |2^, Q are 
discussed in Section ^ Formally of higher order in a, they are numerically 
not negligible. A numerical analysis of the widths and the energy-level shifts is 
carried out in Section |^ at ^^(p*, e^p^) in the chiral expansion. 



2 General features of tttt and ttK atoms 

In this section, we describe the general features of the systems that we are going 
to study. The tt+tt" and iv^K^ atoms are highly non-relativistic, loosely bound 
systems, mainly formed by the Coulomb interaction. The average momentum 
of the constituents in the cm. frame lies in the MeV range. Further, their decay 
widths ~ 0.2eV are much smaller than the binding energies ~ lO'^eV involved. 
The 7r"'"7r~ atom in its ground-state decays predominantly into a pair of two 
neutral pions, through the strong transition 7r"'"7r~ tt'^tt'^. The decay width 
into two photons is suppressed by the factor 4 • lO^'^ ^ For a detailed 
discussion of the decay channels of pionium, we refer to [ p3| . The decays of 
the Tr~ atom have to conserve strangeness. Apart from the dominant S- 
wave decay channel into tt'^X", the only allowed decays are therefore K'^ + 
ni7 -|- n2e~'"e~ and tt'^K'^ + nij + n2e^e~ , where ni -\- n2 > 0. Here ni and n2 
denotes the number of photons and e^e~ pairs, respectively. In the relativistic 
theory, the odd intrinsic parity process + 27 corresponds to a 

local interaction in the Wess-Zumino-Witten term while the transition 

7r"A'+ A:° -f 7 occurs not until 0{p'^ ) H- 

The non-relativistic framework |2^, ^ |24 |25|] we are going to apply, provides 
a systematic expansion in the isospin breaking parameter 5. In the case of 
pionium, we count a as well as (?7i„ — mdY as small quantities of order 5. As 
for the n^K^ atom, both a and rriu — md count as order 5. The different power 
counting schemes are due to the fact that in QCD, the chiral expansion of the 
pion mass squared difference = ~ -^^0 is of second order in to„ — m^, 
while Aif = M"^^ — M^o is linear in m„ — m^. At leading and next-to- leading 
order in isospin symmetry breaking, the tt+tt", (tt^AT"'") atom decays into tt'^tt'^ 
(tt^A'") exclusively. The leading order term for the width is of C((5^/^), isospin 
breaking corrections contribute at order 5^/"^ . The results for the S-wave decay 



widths at order 5^/'^ are presented in Sees. 4.1 and 4.3. At order J^, also other 



decay channels contribute. In Section 5.2, we estimate the order of the various 
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decays. 

The energy-level splittings of the tt+tt^ and tt^K^ atoms are induced by both 
electromagnetic and strong interactions. At order 6^, the energy shift con- 
tributions are exclusively due to strong interactions, while at order S"^, both 
electromagnetic and strong interactions contribute. It is both conventional and 
convenient to split the energy shifts into a strong and an electromagnetic part, 
according to[| 

AE^i = AEl^^ + A£X. (2.1) 

The expressions for the strong energy shift AE^^ at next-to-leading order in 
isospin symmetry breaking are presented in Sees. 4.S and 4.3. The electro- 



magnetic part AE'^ is discussed in Section 5.1. Another important correction 



is generated by the vacuum polarization of the electron. Formally, the vacuum 
polarization contributes to the energy shift at order and to the width at order 
(5^^/^, but these corrections are amplified by powers of the ratio ii+fnie- Here 
/i+ denotes the reduced mass of the bound system and rUe the electron mass. 
The vacuum polarization contributions are discussed in Section |^. 
In what follows, we proceed systematically and discuss in detail the decays and 
bound state spectra within the non-relativistic framework. 



3 Non-relativistic framework 
3.1 Hamiltonian 

The Hamiltonian consists of a infinite series of operators with increasing mass 
dimensions - all operators allowed by gauge invariance, space rotation, parity 
and time reversal must be included. However, in the evaluation of the decay 
width and the strong energy shift at next-to-leading order in isospin symmetry 
breaking, only a few low dimensional operators do actually contribute. For the 
TT^K^ atom, the following Hamiltonian achieves the goal: 

H = Hq + He + Ho + Hs , 

Hr = J d^xnr{0,x), r = 0,C,D,S, (3.1) 

with 

Z— it .0 ^^TT ^ ^ ^ 

= -Cin^_Kln_K^ - C2 {niKlTT^Ko + h.c) - CanlKlw^K, 
-Ci {iriAKlTToKa + h.c) - C5 {TriAKlTr_K^ + h.c) 
-Ce{niTT_)AiKlK^) - C7 {VTT^KlVnoKo + h.c) +■■■ , (3.2) 

-"^Note that this sphtting cannot be understood literally, i.e. there are contributions from 
strong interactions to AE'^. 
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where wAw = uAv + vAu. We work in the cm. system and thus omit terms 
proportional to the cm. momentum. The basis of operators with two space 
derivatives is chosen such that none of them contributes to the S-wave decay 
width and energy shift at the accuracy we are considering. For this reason, we 
transformed the operator with two space derivatives on the neutral fields by the 
use of the equations of motion, 



(3.3) 



For the moment, we further neglect transverse photon contributions. To the 
accuracy we are working, they do not contribute to the decay width and to 
the strong energy-level shifts. However, transverse photons do contribute when 
we work out the electromagnetic energy- level shifts in Section 5.1. The non- 
relativistic Lagrangian in the presence of transverse photons is given in appendix 
The Hamiltonian in Eq. (3.1) is built from the non-relativistic pion and kaon 
fields 



7ri(0,x) 



y di.(p)e'P^a,(p), X,(0,x) -y"di.(p)e'P^b,(p), z = ±,0. 



with dv{i>) ^ (i^p/(27r)^ and 

ai(p),a[(p') 
b,(p),bl(p')' 



= (2^)3<53(p-p')5,fc, 



The two particle states of zero total charge are defined by 

|pi,P2>+ = aL(pi)bt_(p2) |0), |p3,P4)o = aj(p3)bj(p4) |0), 
and the total and reduced masses and /^^ respectively, read 



Mi 



i = +,0. 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



3.2 Master equation 



To evaluate the decay width and the strong energy shifts we make use of resol- 
vents. This technique, which was developed by Feshbach a long time ago p7| , 
has been discussed extensively in Ref. psf . To remove the center of mass mo- 
mentum from the matrix elements of any operators R, we introduce the notation 

a(pi,P2 I R(2) |P3,P4)b 

= (27r)353(p^+p2_p3_p4)Jp^^p2| R(^) |P3^P4)^^ (3.8) 

where a, b stand for 0, +. Further, we have 

a(q, -q I R(z) I p, -p)6 - a(q | B.{z) | p)b. (3.9) 
The master formula to be solved is given by the following eigenvalue equation, 

z-E^- j dv{V){^.nO. P I t{z) I ^'„o, 0) = 0, (3.10) 
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where En — — a^/i+/(2n^) denotes the n-th Coulomb energy and the un- 
perturbed rt-th eigenstate is given by 

I *„o, P> = 1 dKq)*no(q) I il^P + q, i^P - q) + . (3.11) 

Here 4'no(q) stands for the Coulomb wave function of the bound tt^K^ system 
in momentum space. The operator r, defined through 

T = V + VG2-f, V = Ho + Hs, (3.12) 

is regular in the vicinity of £"„. The quantity G" stands for the n-th energy 
eigenstate singularity removed Coulomb resolvent, 

G2 = Gc|l-ydi^(P) |^'„o,P)(*„o,P| |, Go ^ ^ ^ . (3.13) 

The master equation presents a compact form of the Rayleigh-Schrodinger per- 
turbation theory. If we insert r iteratively into ( 3.10| ), the eigenvalue equation 
becomes 

z^En- |*„o(x = 0)|2 [Ci + ClMz)] +■■■ (3.14) 
where ^„o(x — 0) stands for the Coulomb wave function in coordinate space 



and Jq denotes the loop integral in Eq. (D.3). The hmction Jq is analytic in 
the complex z plane, except for a cut on the real axis starting at z = Sq. The 
imaginary part of Jq has the same sign as Im z throughout the cut z plane, 
which does not allow Eq. ( |3.14|) to have a solution on the first Riemann sheet. 
However, if we analytically continue Jq from the upper rim of the cut to the 
second Riemann sheet, we find a solution at z — Rez -I- ilmz, with 

Rez = En + lmz = - ^ , ' Ci^+-- - , (3.15) 



and pn = 2/io {En — So)- In the following, we evaluate the S-wave decay width 
r„o = — 2Imz at order 5^^^ and the energy shift Ai?„o = Rez — £"„ at order 
5^. We focus on the strong part of the energy shift only, the evaluation of the 



electromagnetic energy shift is discussed in Section 5.1. As in Ref. [E3|, we 



reduce Eq. (3.1C) to a one-channel problem with an effective potential W, 

grg = gWg+ gWgG^QTg, (3.16) 

where 

W = V + VpoG" {1 - goygo&cV goV. (3.17) 
Here g, (go) denotes the charged (neutral) two-particle projector 

£' = yrf'^(Pi)rf'^(P2) |pi,P2) + +(pi,P2|, gQ = l-g. (3.18) 
The matrix element of W between charged states takes the formQ: 
+(q| W(z) |p)+ = (27^)3<53(q_p^ \ P' f ^ 



(3.19) 



+w{z) + wi(z)p^ -t- W2{z)q^ + W3(z)pq + 



^The delta function term contributes to the electromagnetic energy shift, see Eq. ( ^.6| ). 
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To the accuracy we are working, only the constant term w{z) contributes to the 
decay width and strong energy shift. We get for the S-wave decay width of the 
TT~K~^ atom at order 5^/^, 

r„o = -2 I *„o(x - 0) |2 lniw{E^a) (1 + 2Rew{E„){SciEn))) + 0{5^), (3.20) 

while the S-wave energy-level shifts due to strong interactions read at order 5^, 

ASj;o = l*no(x-0)|'Rezi;(£;„)(l + Reu;(£;„)(g^'(^„))) + 0(,5^). (3.21) 

The quantity (g" (£'„)), given in appendix is related to the integrated Schwin- 
ger Green function WSl. The real and imaginary part of w{z) are given by 



Rew{En) 
lmw{En) 



-^1 + 4^2 P" 



27r 2 



Cinlpn 5^0 Pn Mi^o 



47r2 



8 M^oMlo 



(3.22) 



The decay width ( 3.20 ) and energy shift ( 3.21 ) still depend on the effective 
couplings Ci, which have to be related to physical quantities. 



3.3 Matching procedure 

We now determine the diverse couplings from matching the non-relativistic and 
t he r ela tivis tic am plitu des at threshold. With the effective Lagrangian in Eqs. 



( |A.2|) , (A. 4) and (|A.5[ ) we may calculate the non-relativistic — > n^K'^ 
and Tr~K~^ Tr~ K'^ scattering amplitudes at threshold at order S. Again, we 
may omit contributions from transverse photons. The radiative corrections to 
the one-particle irreducible amplitudes, generated by transverse photons, vanish 



The coupling C3 enters the decay width ( 3.2C ) at order S^^^ and the energy 



shift ( 3.21 ) at order 6'^ and is therefore needed at 0{S'^) only. However, we 
have to determine both C'l and C2 at next-to-leading order in isospin symmetry 
breaking. The relativistic amplitudes are related to the non-relativistic ones 
through 

r^^''(q;p) -4[a;,(p)c.,(p)c.Kq)^™(q)]^T;ir'(q;P), (3-23) 

where uJi{p) ~ {Mf +p'^y^^ and p, (q) denotes the incoming (outgoing) relative 
3- momentum. In the isospin symmetry limit, only the lowest order of the non- 
relativistic Lagrangian contributes at threshold and the effective couplings Ci, 
C2 and C3 yield 

= — (a+ + a-), C2=-^a^, C, = —a+, (3.24) 

where and Oq denote the isospin even and odd S-wave scattering lengths, 
the notation used is specified in appendix ^ The ttK scattering lengths are 
defined in QCD, at ni^ = rud = m and A/tt = -^tt+j = Mj^+. 
To match the coupling C2 including isospin symmetry breaking effects, we cal- 
culate the real part of the non-relativistic tt^K^ — > tt^K'^ scattering matrix 
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Figure 1: One-photon exchange diagrams for the tt^ — > tt^K'^ scattering 
aniphtudc. The dotted hne denotes a Coulomb photon. The dots stand for the 
coupHngs Ci, i = 1,2. 



element in the vicinity of the threshold at order 6. In absence of virtual pho- 
tons, the real part of the amplitude at threshold reads 

RerOO^±(q; p) = + CaCf Jo(S+)2 + • • • (3.25) 

The ellipsis denotes terms which vanish a threshold or are of higher order in the 
parameter 6. The one-loop integral Jq is given in appendix ^ Bubbles with 
mass insertions and/or derivative couplings do not contribute at threshold at 
order S, since they contain additional factors of and/or S+ — So- 
We now include the Coulomb interaction. Feynman graphs, with a Coulomb 
photon attached such that the heavy fields must propagate in time to connect 
the two vertices, all vanish. This is because we may close the integration contour 
over the zero-component of the loop momentum in the half-plane where there 
is no singularity in the propagators. One example is the self-energy diagram, 
which vanishes at order a. As a result of this, there is no mass renormalization 
in the non-relativistic theory and the mass parameters M^^i and M^., i = 0,+ 



in the non-relativistic Lagrangian (A.4) stand for the physical meson masses. 



The amplitude at threshold contains both infrared and ultraviolet singularities, 
coming from the one-Coulomb photon exchange diagrams depicted in Figure |l]. 
Around threshold, we get for the one- Coulomb exchange diagrams, 

7;°^^^(q; P) = -C2Vcip, P°hr) [1 + ClJ+{P^ij] + C^C2Ba{Pl,) + • • • , (3.26) 

where P^^n: — ^+ + P^/(2m+)- The Coulomb vertex function Vc in Figure |l|(a) 
and the two-loop integral B^, in Fig. ^(c) are given in appendix ^ The integral 



J+, specified in Eq. (D.4), has to be evaluated at d ^ 3, because the vertex 
diagram generates an infrared singular Coulomb phase |^ at threshold. We 
split of this phase ^c, according to 

T;°*(q.p) = e--.T»±(q,p), 

where /i denotes the running scale. The remainder r^R is fr^e of infrared 
singularities at threshold, at order 5. We find for the real part: 

Ref°°^±(q,p) = § +B2lnM + ^ Re<;± +0(p), (3.28) 
IPI M+ 



with 



2 

B^^C2^^ + o{5), B2^~CiC2^ + o{d), (3.29) 
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Figure 2: Non-relativistic ir^ scattering amplitude. The blob de- 
scribes the vector form factor of the pion and kaon. T",^'^ denotes the truncated 
amplitude. 



and 



(A 



K 



At order 6, the constant term in Eq. (3.28) reads 

,3 



= Co 



27r2 



Ci- 



27r 



A. 



1 - A(/x) - In 



(3.30) 



4mI 



(3.31) 



The ultraviolet divergence A(/i), given in Eq. ( |C.5| ), stems from the two-loop 
diagram Be and may be absorbed in the renormalization of the coupling C2, 



1 



2tt 



CiA(^) 



(3.32) 



We now determine the coupling constant Ci. At a = 0, the real part of the 
non-relativistic t:^ scattering amplitude reads at threshold, 



ReT±„'±(p, p)=Ci+ C3C2Vo(S^ 



(3.33) 



where the ellipsis denotes contributions which vanish at threshold or are of 
o{S). In the presence of virtual photons, we first have to subtract the one- 
photon exchange diagram from the full amplitude, as displayed in Fig. ^j. The 
coupling constant Ci is now determined by the one-particle irreducible part of 
the amplitude. The truncated part Tj^^ again contains one-photon exchange 
diagrams as shown in Figure ^, 



f^^iP, P) - -2CiFc(p, PL) [1 + CiJ+iP^)] + C^BaiPL 



(3.34) 



All diagrams with a Coulomb photon exchange between an incoming and an 
outgoing particle vanish, because the pions (kaons) must propagate in time in 
order to connect the two vertices. Again the vertex function Vc leads to an 
infrared singular Coulomb phase at threshold. 



p,p)=e2-«=r±,±(p,p). 



±;±f 



(3.35) 



where T^r^ is free of infrared singularities at threshold at order S. Further, the 



real part of the infrared regular amplitude T^jl^ is given by 



Rer±^±(p,p) 



IpI 



1 



(3.36) 
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Figure 3: One-photon exchange diagrams for the truncated tt^K^ 
scattering ampUtude. The dotted hne denotes a Coulomb photon, 
denotes the couphng C2. 



TT~K+ 

The dot 



with 



and 



B[=CiaTrfi+ + o{S), B'^ = + o(J), 



(3.37) 



-Re A: 



±;± 



= C 



4^^ 



(3.38) 



Here, the ultraviolet pole term A(/i) in Be is removed by renormalizing the 
coupling Ci, according to 



c;(p) = c, 



(3.39) 



The above renormalization of the low-energy couplings Ci and C2 eliminates at 
the same time the ultraviolet divergences contained in the expressions for the 
decay width ( [3.2(]| ) and the energy shift (3.21). We assume that the relativis- 
tic ttK amplitudes at order 5 have the same singularity structure as the non- 
relativistic amplitudes and use Eq. ( 3.23| ) to match the non-relativistic expres- 
sions to the relativistic ones. The calculations of the relativistic tt^ n^K^ 
and — > 'k~ scattering amplitudes have been performed at 0(p'^ , e^p^) 
in Refs. ^ |5^. Both the Coulomb phase and the singular term '-^ ln|p| 
are absent in the real part of the amplitudes at this order of accuracy, they 
first occur at order e^p'*. The quantity ReA™;."*^, (ReA^'^i'') is determined by 
the constant term occurring in the threshold expansion of the corresponding 
relativistic amplitude. Further, the relativistic calculations con- 
tain the same singular contribution ~l/|p| as the non-relativistic amplitude in 
Eqs. (1^ and ( ^^ . 

The results for the matching of the coupling constants C[ (/i) and C2 (/i) yield 
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at next-to-leading order in isospin symmetry breaking, 

2(E+-Eo)(aJ)' 



In 



4^ 



- 1 (a+ + ) 



and 



1 



ReAf^^ + 87r(E+ - So)aJ(ao )' 



-27rQ! In 



-\2 



1 (a+ + ag ) 



(3.40) 



(3.41) 



4 Strong energy shift and width 



The matching relations in Section 3.3, allow us to specify the results for the 
decay width and the strong energy shift in terms of the relativistic t:K scattering 
amplitudes at threshold. The expressions are valid at next-to-leading order in 
isospin symmetry breaking, and to all orders in the chiral expansion. 



4.1 S-wave decay width of the nK atom 



The matching results in Eqs. (3.24) and ( 3.40| ) yield for the decay width at 
order S^^^ in terms of the relativistic Tr^if + n'^K^ amplitude at threshold, 



^ = ^^Re<f + 



o(^). 



(4.1) 



where 



-4:an+{a^ + 



^(n)-V(l) 



1 ^ a 
- + In- 
n n 



+ o{S). (4.2) 



and il;{n) = r'(n)/r(n). Aside from the kinematical factor p*, the decay width 
is expanded in powers of a and m„ — m^;. The outgoing relative 3-momentum 



1 



Pn 



2En 



\{ElMlo,Ml, 



,1/2 



(4.3) 



with A(x, y, z) = x + y + z — 2xy — 2xz — 2yz, is chosen such that the total 
final state energy corresponds to the n-th energy eigenvalue of the tt~ atom. 
In the isospin limit, the Tr~ tt'^K'^ amplitude at threshold is determined 
by the isospin odd scattering length Uq 



In order to extract Oq from the above 
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result of the decay width, we first have to subtract the isospin breaking contri- 
bution from the amplitude. We expand the normalized amplitude in powers of 
the isospin breaking parameter S, 



A — ar, 



0(6). 



(4.4) 



The isospin breaking corrections e have been evaluated at 0{p'^,e'^p'^) in Refs. 
[ pO| pT[ . See also the comments in Section ^ We may now rewrite the expression 
for the width in the following form: 



nO 



2e 



(4.5) 



The corrections 5K,n to the Deser-type formula have been worked out at order 
a, m„ — TOd, arh and {rriu — md)rh. The corrections 5^,„, where n = 1,2 are 
given numerically in Table |[ 

For the P-wave decay width into tt^K^ , the leading order term is proportional 
the square of the coupling C7 and of order 5^^/'^. After performing the matching, 
we get at leading order 



(n' - 1) 



M+P (Ol ) , 



(4.6) 



where denotes the P-wave scattering length. 



4.2 Strong energy shift of the ttK atom 

With the matching results in Eqs. (3.24) and ( ^.41 ), we may specify the S-wave 
energy shifts at order (5**, in terms of the relativistic truncated t:~ K'^ tt' 
amplitude at threshold. 



(4.7) 



with 



1 a 

ipin) - ?A(1) h In- 

n n 



0(6). 



(4.8) 



For the ground-state, the result agrees with the one obtained for the strong 
energy shift in pionic hydrogen p5[, if we replace fi+ with the reduced mass 



of the TT~p atom and ReA^^^^ with the regular part of the n~p amplitude at 
threshold. 

In the isospin limit, the normalized amplitude A' reduces to the sum of the 



isospin even and odd scattering lengths 
powers of a and to„ — to^. 



A' 



+ Oq + e' + o{S). 



Again, we expand A' in 



(4.9) 



The corrections e' have been evaluated at 0{p'^, e^p^) in Refs. See also 

the comments in Section |^. The isospin breaking corrections to the Deser-type 
formula read at next-to-leading order: 



2a V+ 



- + (4.10) 



12 



where E'^ „ has been worked out at order a, m„ — m^, am and (to„ — md)rh m 
the chiral expansion. For the first two energy- levels, the numerical values for 
0^, „ are given m Table |. What concerns the energy splittings for I — 1 and 
n > 2, the leading order contribution is given by 

AEll, = -2C6V*:i(x - 0)V'J'„i(x = 0). (4.11) 

Here denotes the Coulomb wave function with angular momentum ^ = 1. 
The low energy coupling constant Cq is determined through the I = 1 partial 
wave contribution to the relativistic tt~ —^ Tr~ K'^ scattering amplitude and 
we find for the energy shift, 

AE^,, = (4 + «r) . (4.12) 

The result is proportional to the combination + of P-wave scattering 
lengths and suppressed by a factor of . 



4.3 Pionium 



For pionium, we adopt the convention used in Refs. p3| , |30| and count a and 
(rriu — rad)^ as small isospin breaking parameters of order 5. The decay width 
and energy shifts of the tt+tt^ atom can be obtained from the formulae (3.20), 
( p. 21 ) and (4.11) through the following substitutions of the masses M 
, M^o A/ttO and coupling constants^, 



Cx ^ Ci, C2 



2(c2 - 2c4A^), C3 — > 2c3, Cq cq. 



(4.13) 



The factor 2 in substituting C3 comes from the different normalization of the 
7r°7r° state |p3,P4)o — ao(P3)'^o(P4)|0)- coupling constant C2, the sub- 

stitution is non-trivial because our basis of operators with two space derivatives 
differs from the one used in Refs. 0. See also the comment in Section 3.1 
The result for the S-wave decay width of pionium reads at order J^/^ , 



9n3 



a^pl^^Al{l + K^^n) , A^=al-al + e^ + o{5), (4.14) 



where 



^0 



1 , a 

h In- 

n n 



1/2 



oiS), 

(4.15) 
ill. The 



and K — M^^/M^o — 1. The quantity is defined as in Refs. [g3| [ 

isospin symmetry breaking contributions have been evaluated at 0(jr, e^p^) 
in Refs. |30[ The corrections e^r are of the order of a and am. This 
is due to the fact that in the tt+tt^ tt^tt" scattering amplitude at threshold, 
the quark mass difference shows up at order (ruu — md)^m only. For the decay 

^The Ci are the low-energy constants occurring in Refs. [sot . 
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width of the ground-state at order 5^/"^, we reproduce the result obtained in 
Refs. 



20 , E3L BOl . Again we may rewrite the formula for the width: 



2^3 
9n3 



Pl.M-4Y{l + 5^,^) + 0{5'), 5^ 



2e, 



(4.16) 



The parameter i5,r,n contains the isospin breaking corrections to the Deser-type 
formula at next-to- leading order. The numerical values for (Jjr.n, with n = 1,2 
are listed in Table |l|. The decay width of the P-states into a pair of two neutral 
pions is forbidden by C invariance. 

The strong energy shift of the tt+tt^ atom at order 6^ yields 



A' 







i>{n) 



, , X 1 , oi 

V'(l) - - + In- 
n n 



+ 0(6), (4.17) 



where is defined analogously to the quantity A' discussed in Section 4.2. The 
isospin symmetry breaking contributions have been calculated at 0{p^, e^p^) 
in Refs. |5^. Again the corrections are of order a and am. At order S^, 
the Deser-type formula is changed by isospin breaking corrections, according to 



^<no = -^^(2a° + al) (1 + S'^J + 0(6'), <5;,„ = + i^;,„. 

(4.18) 

For the first two energy-levels, the numerical values for S'^ „ are given in Table 
|l]. Finally, the leading order contribution to the strong energy-level shift for 

I = 1 and 71 > 2 reads 

Ai?^,„, = J-^!^a'A4,al (4.19) 
here a\ denotes the P-wave tttt scattering length. 



5 Transverse photons 

We now concentrate on the contributions coming from transverse photons. At 
order S'^, the energy-level shifts in tt+tt" and tt^K^ atoms contain apart from 
the strong energy shift also an electromagnetic contribution as well as finite 
size effects due to the electromagnetic form-factors of the pion and kaon. We 
further discuss the contributions from transverse photons to the decay width 
of the 7r~K^ atom and show that they do not contribute at order S^^^. For 
pionium, the various higher order decay channels have been discussed in Ref. 




5.1 Electromagnetic energy-level shifts 

As mentioned in Section |2| we split the total energy shift AEni in Eq. (2.1) 
into the strong part displayed in Eq. (4.7) and an electromagnetic contribution 
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Figure 4: Self-energy (a) and one-photon exchange (b) to the electromagnetic 
energy shift. The twisted lines denote transverse photons. 



AiS^™. The electromagnetic part is of order and contains both pure QED 
corrections as well as finite size effects due to the charge radii of the pion and 
kaon, see appendix |^. The energy shift contributions due to pure QED have 
been evaluated by the use of the Bethe-Salpeter equation |Q, the quasipotential 
approach ^ and an improved Coulomb potential |Q. Nevertheless, we find 
it useful to provide the calculation within the non-relativistic framework. 
Again, we start with the master equation ( 3.10| ), but instead of the effective 
potential W, we consider the operator f in the second iterative approximation. 



-f = V + VG'^5V + C'(V3). 



(5.1) 



The non-relativistic Lagrangian including transverse photons (A. 2), (A. 4) and 
( A. 5 ) gives rise to the following perturbation, 



V = Hr 



Til 



iA 



1 



ttI V7r_ 



1 



M 



K + 



The photon field A is given by 



(5.2) 



A(0,x) 



(27r)32fcO 



A=l,2 



(5.3) 

where fc" = |k| and e(fc, A) denote the transversal polarization vectors. The 
operator satisfies the commutation relation, 

[a^(fc,A),at(fc',A')] ^ 2k'>{27rfS'ik - k')Sxx' , (5.4) 

and the one-photon states read 

\k,X)^a\{k,X)\0). (5.5) 

The electromagnetic contributions to the energy-level shifts consists of 



1 



1 - 



J rf3r^:,(r)A2*„,(r) + e^X\^no{^ = 0)|^ 



^2 



dj.(p)dz.(p')4':i(p)G^(p,p')^ni(p'), 



(5.6) 
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with fin = + p^/(2/^+) — En and ^ni denotes the Coulomb wave function 
for arbitrary n and /. Here, the first term contains the mass insertions H^, the 
second describes the finite size effects due to Ha , while the last two terms come 
from the self-energy and one-photon exchange diagrams depicted in Figure ^ 
(a) and (b). 

To avoid contributions from hard photon momenta, we use the threshold expan- 
sion |5^, |56| to evaluate the self-energy contributions. This procedure is outlined 
in appendix ^ a nd th e threshold expanded self-energy E/i, where h = Tr,K, is 
specified in Eq. ( |D.17| ). As can be read off from the wave function in momentum 
space, the relative 3- momentum p is of order S. Hence the quantities T,Tr{fin, p) 
and Ejc(^^n, p) count as order and are beyond the accuracy of our calculation. 
What remains to be calculated is the one-photon exchange contribution]^. The 
integrand 



G^(p,p') = 



1 



P-P'l 
(P 



1 

^4 



(p2- 



/2\2 



(p-pO' 



Ip-p 



(5.7) 



is a inhomogeneous function in the parameter 6, and to the accuracy we are 
working required at leading order in S only. 



G^(p,p') = 



1 



1 



4|p-p 



l\2 



(P + P'f- 



2n2 



(P - P')' 



(5.8) 



In order to evaluate the one-photon exchange contributions, we replace the terms 
~ p2^*;(p) and ~ p'^4'„;(p') by making use of the Schrodinger equation. 



r*n;(q). 



(5.9) 



Further, we use the Fourier transform of |p — p'|^^, |p — q|~^ and |p' — ql^^ 
to express the wave functions in coordinate space. The one-photon exchange 
contribution now reads at order a'*, 



|*no(x = 0)|2 



where the expectation values are defined as 



3a^ 
2S^ 



(r- 



/fc = 1,2. 



(5.10) 



(5.11) 



The electromagnetic energy shift at order yields 



4 



+ 



A 9 r 



3/i+ 



_3_ 

8n 



1 



2/ + 1 
3 



1 



7i3(2; + l) 



O(a^lna). (5.12) 



*We thank A. Rusetsky for a very useful communication concerning technical aspects of 
the calculation. 



16 



Figure 5: Annihilation diagram tt+tt s- tt+tt at order a in the relativistic 
theory. 



Here, the first terms is generated by the mass insertions, the second contains 
the finite size efi'ects and the last stems from the one-photon exchange contri- 
butions ( 5.10| ). For n and / arbitrary, we get the same result for the pure QED 
contributions as Ref. |43|. Further the formula for the ground-state 

agrees with the result obtained in Ref. for the electromagnetic energy shift 
of the Tr~p atom, if we replace A through the corresponding quantity in pionic 
hydrogen. 

To analyze the electromagnetic energy splittings of pionium, we need to con- 
struct an effective Lagrangian that describes the relativistic tt+tt^ — > tt+tt^ 
amplitude at threshold correctly up to and including order a. The annihilation 
graph showed in Figure |^ corresponds to a local four pion interaction in the 
non-relativistic Lagrangian. However, the corresponding relativistic scattering 
matrix element vanishes at thres hold. We may thus obtain the electromag- 
netic energy-level shift from Eq. ( ^.12 ) by simply substituting /i+ — > M^+/2, 
S+^2M,+ andA^l/3(r2^), 



Sio 
8n^ 



1 



11 



64^4 2n^{2l + l) 



3(0 



(5.13) 



5.2 Decay channels of the ttK atom 

Next we discuss the contributions from other decay channels to the decay width 
of the 7r~iir+ atom. As already mentioned in Sectional for S-states the only 
possible decay channels are + ni7 + n2e+e^ and ttK'^ + nij + n2e^e~ , 
where ni + n2 > 0- Here ni denotes the number of photons and n2 the number 
of e~^e~ pairs. The decay widths into t:^ + ni7 -I- n2e~^e~ vanish due to 
lack of phase space. Moreover, radiative transitions^ with the emission of one 
photon are forbidden between two states with / = 0. For the 2P-state of the 
TT^ atom on the contrary, the main annihilation mechanism is the 2p — Is 
radiative transition into the ground-state, followed by the decay into tt'^K'^. 
To investigate the decays int o -I- nj, n = 1, 2 we have to extend the La- 
grangian in Eqs. (A^) and (A^) to include terms with odd intrinsic parity, 
such as 



Ca = eDiB ■ {tt[I)KIKo + h.c) -I- e^L'aE • B(7r1; A'^i^o + h.c) 



••, (5.14) 

^Transitions between S-states with the simultaneous emission of two photons are not for- 
bidden. However, they are suppressed by a factor of a^. 
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where uDv = uDv — vDu. The covariant derivative D is specified in appendix 
E denotes the electric and B the magnetic field. The couplings Di and D2 
are real and may be determined through matching with the chiral expansion 
of the relativistic amplitudes. In the relativistic theory, the tt^K^ vertex 
is contained in the Wess-Zumino-Witten term |4^. The t:'^ K'^j interaction 
occurs not until the odd intrinsic parity sector of the ChPT Lagrangian at 0{p^) 
. The such extended Hamiltonian is hermitian and the operator r obeys the 
unitarity condition, 

r(z) - f^z) = -27rif {z)d{z - Hq - lic)T^z). (5.15) 

The symbol 6 is understood as follows: in order to evaluate the right-hand side 
of the equation, we insert a complete set of eigenstates (Hq + Hp) | (3) = Ep | /3) , 
omitting the n-th Coulomb eigenstate of the 'k' atom. This implies for the 
total decay width: 

r„; = ^r„i_/3, (5.16) 



where 



Tni,p = j dpp diy(P)27:diz - Efj){^rd , P | t(z) | /?) 

x(/3|rt(^)|vi/„;,o), (5.17) 



and z is the solution of the master equation (3.10). Here, dpp denotes the 
phase space integral over the intermediate state |/3). At the accuracy we are 
considering, we may use z = £„. In the following, we estimate the order of the 
various decays using this formula. As an illustration, we start with the decay 
into 7r°i^°. The relative 3-momenta of the tt~K~^ pairs p and p' count as order 
d and we have 

dz.(p)di.(p')Ki(p)*„,(p') = OiS^). (5.18) 

As can be read off from the energy delta function, the outgoing tt" and K'^ 3- 
momenta P3 and p4 count as order d^^^ . This leads to a phase space suppression 
factor of order (5^/^, 

diy{p3)diy{pi)d^iP3 + P4)<5 (E^ ~ E^o^o) = 0{5^'^), (5.19) 

where ^ ^ 

E^aK^> ^^0 + + (5.20) 

Further, the reduced matrix element +(p | T{En) \ Pd)o is of 0(1) and the S- 



wave decay width thus starts at order (5^/^. The relation (5.17) allows one to 
rather straightforwardly rederive the next-to-leading order formula for the decay 
width of the S-states. In order not to interrupt the argument, we relegate the 
relevant calculation to appendix ^, and continue here with the discussion of the 
radiative decay into tt^K'^ -\-^. The outgoing tt*' and K'^ 3-momenta again count 
as 0(5^/^), while the outgoing photon 3- momentum k is of order 5. This can 
by seen by performing the phase space integrations over ps, p4 and k explicitly. 
In total, the phase space suppression factor amounts to 5^^"^, 

r/^k 

dv{p:i)dv{pi)-—5''{p:i + P4 -f k)5 (£;„ - E^oko^) = OiS"'^), (5.21) 
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K+ K+ K° 

c d 



Figure 6: Leading order contributions to the decays into 1:^X^-1-^ and + n'^, 
n = 1,2. The dot denotes the couphng C2, while the box stands for the couphngs 
Di and D2 in £a- The twisted Unes denote transverse photons. 



ith 



E. 



Pi 



p1 



(5.22) 



The leading order contribution stems from Figure ||(a) and ^(b). The corre- 
sponding reduced matrix element is given by 

+ (p,-p| VG"V|p3,P4,k,A)o - -eC2p-e(fc,A)[/(M^+,p2,-p-k) 

+ /(A^,+ ,p^p•k)] +••• , (5.23) 



where 



/(M,p2,±p.k) = 1: 



2^2 



_k^ 

2M' 



Ik|± 



p k 



M 



(5.24) 



This matrix element is of order J^/^ which implies that the decay width T-j^oxo^y 
starts at order 6^^/'^. However, this contribution vanishes after performing the 
integrations over p and p'. 

Next, we consider the decay into K'^ + 717, n = 1, 2 (see Figure ^c) and ^(d)). 
Here, the outgoing and photon 3-momenta belong to the hard scale and thus 
count as 0(1). For tt^ A'+ K'^ + 7, the Lagrangian ( ^.14 ) leads to a reduced 
matrix element of order (J"^/^, 



+ (p, -p I V I p4, k, A)o = 2eDip • (k x e{k, A)) 



(5.25) 



Naive power counting implies that the decay width into if" + 7 starts at order 
(5^. The matrix element ( 5.25 ) is odd in p and the S-wave decay width therefore 
even more suppressed, while r„i^o^ starts at order 5^ . 
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For the transition tt K'^ + 27, we get from the Lagrangian ( ^.14 ) a local 

matrix element of order 5, 



(P, -P I V I p4,ki, Ai,k2,A2)o 



e2i:i2fc?e(fci,Ai) • [k2 X e(fc2,A2)] 
-2e^Die{ki,\i) ■ [kj x e(fc2,A2)] 
+fci ^ fc2,Ai ^ A2, (5.26) 



and the decay width of the S-states into + 27 thus starts at order 5^ . For 
the P-wave decay width into + 27 this contributions vanishes, because the 
matrix element in Eq. ( ^.26 ) is p independent. 

Processes with a higher number of photons may be treated in a analogous man- 
ner. We expect them - using power counting arguments - to be even more 
suppressed. Since hard processes such as K*^ + ni7 + ^26+6^, rii + 712 > 
do not contribute to the decay width at order J^^'^, we m ay a ssume tha t all 
couplings in the non-relativistic Lagrangian in Eqs. (A_^), {kA) and (A^) are 
real. The total S-wave decay width of the -k^ K'^ atom amounts to 



r„ 



nO,7roi^o 



(5.27) 



The ■K~ atom in the 2P-state on the other hand decays predominantly 
through the radiative transition into the ground-state. To evaluate this tran- 
sition, we insert the ground-state plus one photon into Eq. ( ^.17 ). Here the 
photon 3- momentum k counts as order 5^, as can be read off from the energy 
delta function 5{E2 — i?i — |k|). At leading order, we get for the spontaneous 
2p — Is transition the well-known expression, see e.g. Ref. |5^ 



r2i = 



(5.28) 



The result is of order and given numerically in Table ^. The first subleading 
decay mode of the 2P-state starts at order 6^ with the odd intrinsic parity decay 
into if ° -I- 7. The P-wave decay width into tt^K^ in Eq. (^) is of order (J^^/^ 
and suppressed with respect to radiative 2p~ Is transition by a factor of 10^''. 



6 Vacuum polarization 



What remains to be added are the contributions coming from the electron vac- 
uum polarization. The calculation of these corrections within a non-relativistic 
Lagrangian approach has been performed in Ref. p2[ . In our framework, the 
contributions due to vacuum polarization arise formally at higher order in a. 
However, they are amplified by powers of the coefficient yu^lrrif.^ where TOe de- 
notes the electron mass. To the accuracy considered here, the only effect of the 
vacuum polarization of the electron is a modification of the Coulomb potential 



Hvac, with 



+ (p| Hvac |q)+ = 5- 



ds 



1 



4m 



2 s -j- (p - q) 



2mi 



1 



4mf 



-,1/2 



(6.1) 



The vacuum polarization leads to an electromagnetic energy shift evaluated in 

AS™^ = (*„n Hvac |*„i). (6.2) 



Refs. 22, E5 M, 
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TT+TT atom 
TT^K^ atom 


5.8 ± 1.2 
4.0 ±2.2 


5.6 ± 1.2 
3.8 ±2.2 


6.2 ± 1.2 
1.7±2.2 


6.1± 1.2 
1.5±2.2 



vac 
nl 



Table 1: Next-to- leading order corrections to the Deser-type formulae. 

For the first two energy-level shifts of pioniun:i|^ and the tt^K^ atom, AE"™ 
and AE^f^ are given numerically in Table || and ^. Formally of order a^'"*"^, 
this contribution is numerically sizeable due to its large coefficient containing 

(M+/me)2'+2. 

The vacuum polarization also interferes with strong interactions and contributes 
to the decay width and to the strong energy shift. This can be seen by insert- 
ing the modified Coulomb potential into the master equation ( 3.10| ). For the 
spectrum and the width of the n^K^ atom, we get 



r„o = ^^PUaof{^ + S.,n + SZ.), 



„3,,2 



^E^no = ^(«f[ + «o)(l+4,„ + '^r„)- (6.3) 

What concerns pionium, the decay width and strong energy shift are modified, 
according to 

A<„„ = -^^(2«^ + «o)(l + '5;,„ + Cn)- (6-4) 

The correction (5™^, h — tt,K is proportional to the change in the Coulomb 
wave function p3| of the bound system due to vacuum polarization, 

2^^„o(x =. 0) 

Here, stands for a generic Coulomb wave function and h = tt, K. For the 
ground-state, this result is contained in Table II of Ref. Formally, (5™^ is 

of order a^, but enhanced because of the large coefficient containing /i-|_/me. 



7 Numerics 

In the numerical evaluation of the widths and energy shifts of the tt+tt" and 
TT^K^ atoms, we use the following numbers: The tttt scattering lengths yield 
aj] = 0.220 ± 0.005, a§ = -0.0444 ± 0.0010 and a\ = (0.379 ± 0.005) • lO'^ Af-+ 
H . The correlation matrix for aj] and ag is given in Ref. Q . For the isospin 



symmetry breaking corrections to the tttt threshold amplitudes (4.14) an d ([4.17 ) 
at order eV, we use = (0.61±0.16)-10-2 and = (0.37±0.08)-10~2 ^^|5j 

®For pionium, the electromagnetic energy shift due to vacuum polarization is denoted by 
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TT^vr atom 










n=l, 1^0 
n=2, 1^0 
n=2, l^l 


-0.065 
-0.012 
-0.004 


-0.942 
-0.111 
-0.004 


-3.8 ±0.1 
-0.47 ±0.01 
~ -1 • 10"^ 


2.9±0.1 
23.3 ±0.7 
~ 1.2 • 10"* 



Table 2: Numerical values for the energy shift and the lifetime of the tt+tt 
atom. 



The values for the ttK scattering lengths are taken from the recent analysis 
of data and Roy-Steiner equations The S-wave scattering lengths yield 

a+ = (0.045 ± 0.012)M-+\ = (0.090 ± 0.005)M-+i|3|, and for the P-waves 

we use a/^ = (0.19 ± 0.01) • lO-^Af-^^ and of ^ = (0.65 ± 0.44) • lO-^M^+^H. 
The correlation parameter for Qq and Qq is also given in Ref. |p3|. The isospin 



breaking corrections to the t:K threshold amplitudes (4^) and ( |4.9| ) have been 
worked out in Rcfs. |5^, pT| at Ojp'^, e^p^). Whereas the analytic expressions 
for e and e' obtained in Refs. |34 |50| , pT] are not identical, the numerical values 
agree within the uncertainties quoted in [5l|] . In the following, we use [pT] 
e = (0.1 ± 0.1) • 10-^ M-^ and e' = (0.1 ± 0.3) • lO^^Af-^i. For the charge 
radii of the pion and kaon, we take {r'^+) — (0.452 ± 0.013) fm^ and {r'^+) — 
(0.363 ± 0.072) fm^ ||5^. In the evaluation of the uncertainties, we take into 
account the correlation between the tttt, (ttK) scattering lengths. 



The isospin breaking correction s Sh,n and /i = tt, X to the widths (4.5), 
(4.16) and strong energy shifts (4.1C), (4.18) are listed in Table |[ The energy 
shift corrections (5^ „ are smaller than in the case of pionium. This discrep- 
ancy comes from the different size of the isospin breaking contributions to the 
elastic one-particle irreducible tttt and ttK amplitudes. At leading order in the 
chiral expansion, the isospin breaking part of the ir^ — > Tr~K~^ amplitude 
at threshold is suppressed by a factor of M^+/AI^+ with respect to the corre- 
sponding quantity in tttt scattering. 

As described in Section ^, Eqs. (3.3) and (6.4), the width and strong energy 



shift are modified due to vacuum polarization, according to 

Sh,n ^ Sh,7i + Sl"^^, <5^,„ (^J,.,,! + (5™^, (7.1) 



where h = tt,K and (5™^ is defined in Eq. (6.5). For the ground-state, the 
corrections due to vacuum polarization yield 6]^\ — 0.45 • 10~^ and (5™^^ = 0.31 • 
10^^||2^. However, for the numerical analysis of the width and the strong energy 
shift, we may neglect the contributions from (5™^, because the uncertainties in 
Sh.n and S'l^ „ are larger than SJ^^^^ itself. 



For the first two states of the tt+tt and n^K^ atoms, the numerical values 
-1 



for the lifetime r„; = ^ni^j {TTT,ni = T^^^) and the energy shifts are listed in 



Table ^ and |[ The numbers for the lifetime and strong energy shifts of the 
S-states are valid at next-to-leading order in isospin symmetry breaking. The 
bulk part in the uncertainties of these quantities is due to the uncertainties 
in the corresponding scattering lengths. For the lifetime of the 2P-state, the 
numerical values are valid at leading order only, and determined by the 2p — Is 



radiative transition in Eq. (5.28) | |59| |. 

The energy-level shift due to vacuum polarization A£;™^ (A£;™^;) g || 
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TT^K^ atom 




Ai?-[eV] 




io^^T,4s] 


n=l, 1=0 
n=2, ;=0 
n=2, 1=1 


-0.095 
-0.019 
-0.006 


-2.56 
-0.29 
-0.02 


-9.0± 1.1 

-l.lzbO.l 

~ -3-10-6 


3.7±0.4 
29.4 ±3.3 
~ 0.7- 10"^ 



Table 3: Numerical values for the energy shift and the lifetime of the n^K^ 
atom. 



is specified in Eq. (6.2). Formally of higher order in a, this contribution is 
numerically sizable. We do not display the error bars for the electromagnetic 
energy shifts. At order a^, they come from the uncertainties in the charge radii 
of the pion and kaon only. In the case of pionium, the uncertainties of AE™Iq at 
order amount to about 0.7%, while for the n^K^ atom AiJ™ is known at the 
5% level. To estimate the order of magnitude of the electromagnetic corrections 
at higher order, we may compare with positronium. Here, the and Ina 
corrections |Q amount to about 2% with respect to the contributions. 
Both, the electromagnetic and vacuum polarization contributions to the energy 
shift are known to a high accuracy. A future precision measurement of the 
energy splitting between the nS and nP states [ |6l| wi ll therefore allow one to 



extract the strong S-wave energy shift in Eq. (4.10), and to determine the 
combination a|j' + Oq of the nK scattering lengths. The 2s — 2p energy splitting 
yields 

Ai?2s-2p = Ai?20 + AiJj™ ^ Ai?2™ + Ai?2o'^ — Ai?2i'^ 

= -1.4±0.1eV. (7.2) 

The uncertainty displayed is the one in AE^ only. To the accuracy we are 
working, we may neglect the strong shift in the 2P state, it is suppressed by the 
power of . For pionium the energy splitting between the 2S and 2P states 
reads 

AE^^2s-2p = A£;^_2o + AS™2o - AE^^^i + A£:™20 - A£;™2i 

= -0.59 ± 0.01 eV. (7.3) 

Again the uncertainty displayed is the one in Ai?Jj 2 only ^-nd we may neglect 
contributions from the strong shift in the 2P state. 

As an illustration, we alternatively use the ChPT predictions for the wK scatter- 
ing lengths HJ, |3^ in the numerical evaluation of the lifetime. The ChPT predic- 
tions yield at order a+ = (0.032 ± 0.016)M-+^ and = (0.079 ± 0.001)Af-/ 
[ p3[ . Here, the errors include the uncertainties in the values of the input param- 
eters only. The uncertainty in Oq is remarkably small, because the isospin odd 
scattering length involves at C'(p^) a single low-energy constant L5 On 
the other hand, contains apart from the combination 2Lg -I- Lg five further 
coupling constants |Q , which are enhanced by one power of Af ^+ /Af^+ with 
respect to the counterterm in Oq . For , the two- loop correction has to be 
rather substantial, in order to reproduce the central value of the Roy-Steiner 
evaluation. Very recently, the chiral expansion of the ttK scattering amplitude 
at next-to- next-to-leading order became available |^ . According to the prelim- 
inary numerical study performed in Ref . [p7[ , the S-wave scattering lengths are 
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at order in reasonable agreement with the Roy-Steiner evaluation Q. The 
0{p'^) ChPT prediction for the lifetime of the tt^K^ atom in its ground-state 
is 

no = 4.7 • 10-'% ChPT, [Oip% (7.4) 

whereas 

no = (3.7 ± 0.4) • 10-^^s, Roy-Steiner. (7.5) 

The ChPT prediction is valid at next-to-leading order in isospin symmetry 
breaking and up to and including 0{p'^) in the chiral expansion. 



8 Summary and Outlook 

We have considered the spectra and decays of tt+tt" and n^K^ atoms in the 
framework of QCD -|- QED. We evaluated the corrections to the Deser-type 
formulae for the width and the energy shift - valid at next-to-leading order in 
isospin symmetry breaking - within a non-relativistic effective field theory. It 
is convenient to introduce a different book keeping for the tt+tt^ and n^K^ 
atoms. What concerns pionium, we count a and (m„ — m^)^ as small quantities 
of order 6, in the case of the tt^K^ atom both a and m„ — rrid are of order 
5. The different counting schemes are due to the fact that in QCD, the pion 
mass difference starts at (to„ — m^Y , while the kaon mass difference is linear in 
TUu - rud- 

Consider first the energy shifts that are split into an electromagnetic and a 
strong part, according to Eq. (2.1). The electromagnetic part in Eqs. ( ^.12 ) 



and ( 5. 13] ) contains both pure QED contributions as well as finite size effects due 



to the charge radii of the pion and kaon. The strong energy shift of the tt~ 
atom is proportional to the one-particle irreducible ir^ — > tt~ scattering 
amplitude at threshold. In the isospin symmetry limit, the elastic threshold 
amplitude reduces to the sum of isospin even and odd scattering lengths UQ+aQ. 
The isospin breaking contributions to the amplitude have bee n eva luated at 
0(e2p^p4)g |Jj tj^j; framework of ChPT. The result in Eq. ( |4.10| ) displays 



the S-wave energy shift in terms of the sum a J -f Cq , and a correction of order a 
and m^—md. For the first two energy- level shifts, the isospin symmetry breaking 
correction modifies the leading order term at the 2% level. The isospin even 
scattering length is sensitive to the combination of low-energy constants 
2Lq + Lg. The consequences of this observation for the SU(3)xSU(3) quark 
condensate py] remain to be worked out. In the case of pionium, the strong 



energy shift displayed in Eq. ( 4.18 ) is related to the tttt scattering lengths 
combination 2ao -I- and a correction of order a and (m„ — iridY ■ For the 
first two energy-levels, these isospin symmetry breaking contributions amount 
to about 6%. 

A future measurement of the energy splitting between the 2S and 2P state in 
the TT+TT", [tt'^K'^) atom will allow one to extract the strong energy shift and 
to determine the scattering lengths combination 2ao-|-ao, {a^ -{-a^). This is due 
to the fact that the electromagnetic energy shifts are known to high accuracy 



and the strong P-wave energy shifts in Eqs. (1.12) and (4.19) are suppressed 
by the power of a^. However, another higher order correction - generated by 
the vacuum polarization - is numerically sizable and contributes to the energy 
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splitting in Eqs. (7.2) and ( |7.3| ). Formally of order a^'+^, this correction is 
enhanced due to its large coefficient containing (/i+/me)^'+^. 
We now turn to the decay widths of the tt+tt" and tt^ atoms. At lead- 
ing and next-to leading order the tt+tt" and n^K^ atoms decay into 27r° and 
TT^K'^ exclusively. Aside from a kinematical factor - the relativistic outgoing 3- 
momentum of the bound system - their decay width can be expanded in powers 
of a and rriu — md- By invoking ChPT, the result for the S-wave decay width of 
the t:~ atom may be expressed in terms of the isospin odd scattering length 
Uq , and an isospin symmetry breaking correction of order a and uiu — md, see 
Eq. (4^). For the ground-state decay width, this correction modifies the lead- 
ing order Deser-type relation at the 4% level. The next- to-lea ding order result 
for the S-wave decay width of pionium is given in Eq. (4.16). The expression 
for the ground-state width agrees with the result obtained in Refs. H^, po|| . 
For the 2P state of the 'k~ , (tt+tt") atom, the decay width starts at ord er 
with the 2p— Is radiative transition into the ground-state, s ee E q. ( ^.28 ). The 
P-wave decay width of the ii^ atom into tt^K'^ in Eq. (4^) is suppressed 
by the power of 5^'^^'^ . For pionium, the P-wave decay width into a pair of two 
neutral pions vanishes due to C invariance. 

We find it very exciting that in view to the beautiful work performed by our 
experimental colleagues, we may expect that many of the above predictions will 
be confronted with experimental data in a not too distant future. This will 
certainly improve our knowledge of the low-energy structure of QCD. 
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A Non-relativistic Lagrangian 

The non-relativistic Lagrangian must be invariant under space rotation, P, T 
and gauge transformations. We do not include terms corresponding to transi- 
tions between sectors with different numbers of heavy fields (pions and kaons). 
These interactions describe processes with an energy release at the hard scale. 
In general, such decay processes are accounted for by introducing complex cou- 



plings in the non-relativistic Lagrangian. However, as shown in Section 5.2, 
intermediate states do not contribute to the decay width at order (5^/^ and we 
may therefore assume the low-energy couplings to be real. 
In the sector with one or two mesons, the non-relativistic effective Lagrangian 
is given by 

= £1 +4"^ +4'^ + ••• (A.l) 
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The first term contains the one-pion and one-kaon sectors 



A A2 N 

V- t / D2 D4 X 

+ > ttI i A - A^Tr+ H 1 ^ H 7r± 

_|_ /i 7j--r 

with electromagnetic charge e, E = — VAq ^ A and B = V x A. The covariant 
derivatives of the charged meson fields are given by 

Dt'K± — dt7T± T ieAQTr±, Dttj- — V7r± ± ieA.'K±, 

DtK± = dtK± =F ieAoK±, T>K± = Viv:± ± ieAK±. (A.3) 

The one-pion-one-kaon sector of total zero charge reads at lowest order 

4") = CnrlKl7r_K^ + C2 {nlKliToKo + h.c) +C37rlKlTT„K„. (A.4) 

To evaluate the decay width and energy shifts of the 7r~7<r+ atom, we need in 
addition the following terms with two covariant space derivatives]^ 

' = C4 (TriQ^KU^K,, + h.c) + C5 (niB^KlTr_K+ + h.c 



+Ce{TTiTT_)B^{KlK\) + C7 (VTrt isrtvTToJs:,, + h.c) + • • • ,(A.5) 

where mD w = mD^w + uD^it. We work in the center of mass system and thus 
omit terms proportional to the cm. momentum. We do not display time deriva- 
tives, for on-shell matrices they can be eliminated by the use of the equations of 
motion. The parameters M^^i (Mj^i) denote the physical pion, (kaon) masses - 
there is no mass renormalization in the non-relativistic theory, see Section 3.3. 
We work in the Coulomb gauge and eliminate the A" component of the photon 
field by the use of the equations of motion. At the accuracy we are considering, 
the term linear in D(E) in Eq. ( [A.2| ) then reduces to a local interaction which 
modifies the low-energy coupling Ci, 

C[ = Ci - e^A, 

7r+ K+ 



^In the cm. system and by the use of the equations of motion, we identify 
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It is sufficient to matcli tlie non-relativistic couplings and Ck at order 5'^. 
We tfierefore consider the pion and kaon electromagnetic form-factors in the 
external field A^. The results of the matching yield 

c^ = Ml^{rl), c,, = Ml^{rl), (A.7) 

where and denote the charge radii of the charged pion and kaon, respec- 
tively. The remaining low energy constants Ci, i — 1 , ... 7 may be determined 
through matching the ttK amplitude in the vicinity of the threshold for different 
channels, see Section b.3l 



B Relativistic scattering amplitudes 

First, we consider the 5=1 processes 

Tr-K+ ^TT-K+, n-K+^Tr°K°, tt^K^^tt^K^, (B.l) 

in the isospin symmetry limit. The decomposition into amplitudes with definite 
isospin yields 

T±'± = \[T'^\s,t) + 2T''\s,t)], 

rO°'± = :/l[T^/^{s.t)-T^'\s,t)], 

T"°'°" = -[2T^'^s.t)+T"''{s,t)]. (B.2) 
3 

The isospin 1 — 1/2 and / = 3/2 components are related via 

T^'\s, t, u) = ^T'^'iu, t, s) ~ ^T'^^is, t, u). (B.3) 

The T+, {T-) amplitude 

T+{s,t) = i[T^/^(s,i) + 2T^/^(s,t)], 

T-{s,t) = i[T^/^(,s,i)-T^/^(s,t)], (B.4) 

is even (odd) under s u crossing. In the s-channel, the decomposition into 
partial waves reads 

oo 

T\s,t) = 167r^(2; + l)if(s)P/(cos6l), (B.5) 
/=o 

where s = [w^+(p) -t-ti;if+(p)]^, t — — 2p^(l — cos6') and 9 is the scattering angle 
in the cm. system. 

The real part of the partial wave amplitudes near threshold is of the form, 

Retf {s) = ^p2' {a\ + 6f p2 + 0(p^)) , (B.6) 
where a\ denote the scattering lengths and h\ the effective ranges. 
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What follows, are the tttt scattering processes: 

7r"7r+ ^ 7r"7r+, 7r"7r+ ^ 7r°7r", (B.7) 

with 

Tf'^ - i[T2(.,t)+3T,i(.,0 + 2r°(.,t)], 
T™^± = \[n{s,t)-n{s,t)\. (B.8) 
The decomposition into partial waves yields 

oo 

Tl{s,t)=?,2^Y.^2l + \)Pi{cosO)ti^i{s), (B.9) 

where s = 4(A/^+ + p^). At threshold, the partial wave amplitudes take the 
form 

Reti^iis) = [a\ + p2&f + 0(p4)] . (B.IO) 



C Schwinger's Green function 

The Schwinger Green function fulfills 



E- 



2^l+ 



+(q| Gc(z) \p)+ + e' / d^fc) 



(k-q) 



^(k|Gc(z)|p)H 



(27r)353(q_p), (c.l) 

where E — z ~ E+. The explicit expression is given by^ 

(27r)3,53(q- p) f 4™ 



,(q|Gc(z)|p)+ = 



E-^ 



1 



s - ^ (q - p)2 s - ^ 



1 



■4TTariI{E, q, p) 



2/i+ 



(C.2) 



with = a/2(-S/(2/i+))-i/2. The function 



J(i?, q,p)= Z'dpp-" 
Jo 



(q-pfp+4(l-p^)(i^-^)(i?-^ 
a'^ 2/i+ z/i+ 



(C.3) 



contains poles at ry = 1, 2, . . . or, equivalently, at z = E^- The integral 



(27r)'* (27r) 



^+(q|GS(i;„)|p)+, 



- ^(1) - - + i 

71 2 



A(/i) - 1 + 21n- + In 



(C.4) 



^To simplify the notation, wc omit the positive imaginary part in E. 
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with ip{n) = T'{n)/T{n) develops an ultraviolet singularity as d — > 3, 

1 



2\d-3 



d-3 



r'(l) - ln47r 



(C.5) 



D Non-relativistic integrals 

What follows is a list of the non- rcla tivistic integrals used to calculate the ttK 
scattering amplitudes in Section 3.3 as well as in the evaluation of the decay 
width (Sec. 3^) and energy shifts (Sees. and jA). Whenever necessary, 
the integrations are worked out in D ^ A space-time dimensions to take care of 
possible ultraviolet or infrared divergences. 
The non-relativistic propagator of the heavy fields reads 



G^HMix) = i{0\Th{x)hHO)\0) = 



1 



(D.l) 



where h stands for the free fields, with h — tt , tt , K+,K" and Mh denotes the 
corresponding mass. The tadpole diagrams C?nr,/i(0) vanishes in dimensional 
regularization. This can be seen, by performing the integration over the zero 
component of the loop momentum explicitly. The remaining integral is scaleless 
and therefore zero in dimensional regularization. 

At a = 0, the elementary loop function to calculate a diagram with any number 
of bubbles is given by 



1 



1 



2M i 



2M„ 



po + io' 



(D.2) 



where i = 0. After the integration over the zero component of the loop 
momentum, we arrive at 



MP") 



1 



Ji. 



(D.3) 



The function is analytic in the complex plane, except for a cut on the real 
axis for P° > S^. For P° > and d 3, we get 



UP") 



3/2 



V27r 

X (—ITT 



d-Z 



r'(l)+ln47r + ln2/i,(pO-S,)) 



(D.4) 



To obtain the contribution to the scattering amplitude, we need to evaluate this 
function at = u)^+{Tp) + ^K+iv)- At threshold, the integral J^{P^) is of 
order p, while Jo(P°) is given by 



(D.5) 
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We now include the Coulomb interaction. The self-energy diagram with one 
virtual Coulomb photon vanishes, because the integration contour over the zero 
momentum of the photon can be closed in the upper half-plane, where there 
is no singularity. Next, we evaluate the Coulomb vertex function Vc and the 
two-loop integral Be- The Coulomb vertex function, see for example Figure 
0(a), is given by 



1 



(27r)^ Ip - ip M.. + ^ 



1° 



2M„ 



(D.6) 



After integrating over the zero component of the loop momentum, the function 
amounts to 



1 



1^ 

2m+ 



(D.7) 



The contribution to the scattering amplitude is obtained for P° — uj^+{p) 
(p) ■ We expand the function around threshold which leads to 



^c(p,Pthr) 



(D., 



where P^^^^ = E+ -I- and the infrared-divergent Coulomb phase Oc is specified 



in Eq. (3.27). The two-loop Coulomb photon exchange diagram depicted in 
Figure ll|(c) reads 



i?c(P") 



d^h d^h 1 

1 1 



1? 



2M^ + 
1 



1? 



1 



- po + ;o 



2M 



1q 



2M, 



(D.9) 



Performing the integrations over the zero components of the loop momenta li 
and I2, we get 



d% d% 1 



1 



1 



{2nr (2^Y I li - I2 P po _ s+ - ^ PO - S+ - ^ ' 

+ 2n+ + 2n+ 



The expansion in the vicinity of the threshold amounts to 

2|p 



Bc{P, 



^ 
tluv 



2tt 



A{fi) + 21n- 



1 ~ in 



0{d~3). 



(D.IO) 



(D.ll) 



The ultraviolet pole term A(/i) is given in Eq. (C.5). 

In the presence of transverse photons, the non-relativistic integrals have to be 
treated properly in order to avoid loop momenta coming from the hard scale. 
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Otherwise, these loop corrections lead to a breakdown of the non-relativistic 
counting scheme. We use the threshold expansion js^, |5^ to disentangle the 
hard scale (given by the meson masses) from the soft scales. We illustrate the 
procedure for the two-point function of the charged pions and kaons at order 



I / dxe'P^O\Th±{x)hi{0)\0) 



1 



(D.12) 



where h — 7t,K and denotes the corresponding mass. The self-energy 'Sh is 
depicted in Figure 0(a) (upper line). In D 4 space-time dimensions, we have 



p2 - (p • k)Vk2 



MliJ {27r)D _k2^M+ + ^ 



fcO) 



+ Oie^). (D.13) 



After integrating over the zero component of the loop momentum, we arrive at 



S„(/,p) 



e 

Ml 



d'^k 1 - (p-k)Vk2 

(2^)'i2|k|^^ 



2Af I 



Pk 

M+ 



(D.14) 



where ft = Af++p^/ (2M+)— p". The threshold expansion amounts to expanding 
the integrand in Eq. ( D.14 ) in the small parameter v, according to the counting^, 



k = 0(w2), p^O{v), n^O{v^). 



The threshold expanded self-energy 



M-2 r(d)r(2-rf) , 4, 



2A/^" (47r)'^/2r(l + I) 

contains an ultraviolet divergence as d — > 3, 



^h[n,p) 



67r2Af2 



p^n L{n) + In 



2n 



L(/i) = 



d-3 



-^4(r'(l)+ln4.+ l) 



(D.15) 
(D.16) 



+ 0(e4,d-3). 



(D.17) 



E Unitarity condition: Evaluation of the width 



The unitarity condition in Eq. (5.17) renders the evaluation of the S-wave decay 
widt h at n ext-to-leading order straightforward. It can be seen fromEqs. (|l8|) 
and ( ^.19| ) that in order to evaluate the width at 0{5^^'^), it suffices to calculate 
the matrix element +(p | f{En) \ P3)o at order 6. Here, the following term occurs 



-c. 



'I 7r" ^ 

2 8 A/3 A/3 „ 



|*no(x = 0)py" di.(p3)27r^(s„ - Eo 



4 


1 


XP3 






z - So - 



Pi 

2^0 



? V P3 



Pj A 

2Aio/ 



(E.l) 



^Instead of first performing the integration over th e zero component of the photon field, 
we may apply the threshold expansion directly to Eq. ( [D.13| ), with fc" = 0{v'^). 
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which is generated by the matrix element +(p | HsGc(2;)Ho | P3)o and its her- 
mitian conjugate. This contribution can be calculated by the use of 



2-Ki5 E, 



V pi 

2^0 

Mo 



1 



1 



z - So 



Pi 

Mo 



So 



Pi 

Mo 



(E.2) 



z^E„+if 



and the result for the decay width at order S^^^ agrees with Eq. (3.20). 
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